1 



Master equations for perturbations of generalised 
static black holes with charge in higher dimensions 

Hideo KODAMA^'*) and Akihiro Ishibashi^'**) 

^ Yukawa Institute for Theoretical Physics, Kyoto University 
Kyoto 606-8502, Japan 
' "^D.A.M.T.P., Centre for Mathematical Sciences, University of Cambridge 

f-*) ■ Wilherforce Road, Cambridge CBS OWA, United Kingdom 

(N : 

< I ■ We extend the formulation for perturbations of maximally symmetric black holes in 

' higher dimensions developed by the present authors in a previous paper to a charged black 

\ hole background whose horizon is described by an Einstein manifold. For charged black 

. holes, perturbations of electromagnetic fields are coupled to the vector and scalar modes 

' of metric perturbations non-trivially. We show that by taking appropriate combinations of 

' gauge-invariant variables for these perturbations, the perturbation equations for the Einstein- 

, Maxwell system are reduced to two decoupled second-order wave equations describing the 

^ ■ behaviour of the electromagnetic mode and the gravitational mode, for any value of the 

QQ ' cosmological constant. These wave equations are transformed into Schrodinger-type ODEs 

^SJ ' through a Fourier transformation with respect to time. Using these equations, we investigate 

( , the stability of generalised black holes with charge. We also give explicit expressions for 

00 ' the source terms of these master equations with application to the emission problem of 

' gravitational waves in mind. 

cn : 
o . 

■ §1. Introduction 

> \ 

pi i' In recent years, motivated by developments in higher-dimensional unification 

1^ ■ theories, the behaviour of gravity in higher dimensions has become one of the major 

subjects in fundamental physics. In particular, the proposals of TeV gravity theories 
' in the context of large extra dimension^^''^ and warped compactificatiorP^''^ have 

^ ■ led to the speculation that higher-dimensional black holes might be produced in 

H ; cohiderPJ^^? ' and in cosmic ray events.®' ^ ' Although fully nonlinear analysis 

of the classical and quantum dynamics of black holes will eventually be required to 
test this possibility by experiments, linear perturbation analysis is expected to give 
valuable information concerning some aspects of the problem, such as the stability 
of black holes, an estimation of gravitational emission during black hole formation, 
and the determination of the greybody factor for quantum evaporation of black 
holes .'^''^'^'^ The linear perturbation theory of black holes can be used also in the 
quasi-normal mode analysiP»'E3JnSinni,in»SlMJJini,EJ,ll3 of the AdS/CFT issues 
and to obtain some insight into whether the uniqueness theorems of asymptotically 
flat regular black holes in four dimension (see Ref. I23|) for a review) and in higher 
dimension^ 'I^J^ 'El -I^J^'EnJ can be extended to the asymptotically de Sitter 
and anti-de Sitter cases. 

With the motivation provided by the developments described above, in a pre- 
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vious paper, Ref. 31) (Paper I), we developed a formulation that reduces the linear 
perturbation analysis of generalised static black holes in higher-dimensional space- 
times with or without a cosmological constant to the study of a single second-order 
ODE of the Schrodinger-type for each type of perturbation. Here, a generalised 
black hole is considered to be a black hole whose horizon geometry is described by 
an Einstein metric. This includes a maximally symmetric black hole, i.e. a black hole 
whose horizon has a spatial section with constant curvature, such as a spherically 
symmetric black hole, as a special case. Then, in Ref. I32p (Paper II), we studied 
the stability of such black holes using this formulation and proved the perturba- 
tive stability of asymptotically fiat static black holes in higher dimensions as well 
as asymptotically de Sitter and anti-de Sitter static black holes in four dimensions. 
We also showed that the other types of maximally symmetric and static black holes 
might be unstable only for scalar-type perturbations. 

One of the main purposes of the present paper is to extend the formulation given 
in Paper I to a generalised black hole with charge and analyse its stability. This 
extension is non-trivial, because perturbations of the metric and the electromagnetic 
field couple in the Einstein-Maxwell system. Hence, the main task is to show that 
the perturbation equations for the Einstein-Maxwell system can be transformed into 
two decoupled equations by an appropriate choice of the perturbation variables, as in 
the four-dimensional case.'^''^''^ Since higher-dimensional unified theories based 
on string/M theories contain various U{1) gauge fields, this extension is expected to 
be useful in studying generic black holes in unified theories. 

The other purpose of the present paper is to give explicit expressions for the 
source terms of the master equations. This information is necessary to apply the 
formulation to the estimation of graviton and photon emissions from mini-black holes 
formed by colliders. 

This paper is organized as follows. In the next section, we first make clear 
the basic assumptions regarding the unperturbed background, and then we give a 
general argument concerning the tensorial decomposition of perturbations. We also 
give basic formulas for the perturbation of electromagnetic fields that are used in 
later sections. Then, in the subsequent three sections, we derive decoupled mas- 
ter equations with a source for the Einstein-Maxwell system in a generalised static 
background with a static electric field for tensor-type, vector-type and scalar-type 
perturbations. In §6, using the formulations given in the previous sections, we anal- 
yse the stability of generalised static black holes with charge. Section 7 is devoted 
to summary and discussion. 

§2. Background Spacetime and Perturbation 

In this section, we first explain the assumptions concerning the unperturbed 
background spacetimes considered in the present paper and present the basic formu- 
las concerning them. Then, we give general arguments on the types of perturbations 
and the expansion of perturbations in harmonic tensors on the Einstein space by 
supplementing the argument of Gibbons and Hartnoll given in Ref. I35|) with some 
fine points associated with scalar and vector perturbations. Finally, we give the basic 
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perturbation equations of electromagnetic fields and formulas used in the subsequent 
sections. 

2.1. Unperturbed background 

In the present paper, we assume that the background manifold has the structure 
M ~ AA^ X /C" locally and its metric is given by 

ds^ = gab{y)dy''dy' + r\y)dal (2-1) 

Here, do"^ = ^ij{z)dz'^dz^ is the metric of the n-dimensional Einstein space /C" with 

Ri^ = {n-l)K^ij, (2-2) 

where Rij is the Ricci tensor of the metric ^ij. When the metric ()2-l() represents 
a black hole spacetime, the space /C" describes the structure of a spatial section 
of its horizon. In the case in which /C" is a constant curvature space, K denotes 
its sectional curvature, while in a generic case, K is just a constant representing 
a local average of the sectional curvature. Because an Einstein space of dimension 
smaller than four always has a constant curvature and is maximally symmetric, this 
difference arises only for n > 4, or, equivalently, when the spacetime dimension 
d = n + 2 is greater than or equal to 6. In the present paper, we assume that /C" is 
complete with respect to the Einstein metric, and we normalise K so that K = 0, ±1. 
The non-vanishing connection coefficients of the metric H2-l[) are 

nc = 'my). r;k = r;dz), n^ = -rD\^i,{z), r;, = ^6), (2-3) 

and the curvature tensors are 

R\cd = 'R\cd, (2-4a) 

R%bj = ^9^J, (2-4b) 

R'jki = R'jki - {Dr)\5l^ji - 6lj,k), (2-4c) 

where R^jki is the curvature tensor of 'jij. From this and H2-2() . we obtain 

l2„ DaDbV 

Rab= 77 Rgab - n , (2-5a) 

2 r 

Rai= 0, (2-5b) 

R,,= i^-— + (n - 1) gi„ (2-5c) 

□r K — (DrV 
R =^R-2n— +n{n-l) — (2-5d) 

where Da, n and ^R are the covariant derivative, the D'Alermbertian and the scalar 
curvature for the metric gab, respectively. Thus, the Ricci tensor takes the same 
form as in the case in which K,^ is maximally symmetric, as first pointed out by 
Birmingham 
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As the background source for the gravitational field, we consider an electromag- 
netic field whose field strength tensor J^^y has the structure 

^ = \Eoeabdy'' A dy^ + ^J^ijdz' A dz^ . (2-6) 
Then, from the Maxwell equation dj^ = 0, we obtain 

Eq = Eo{y), Tij = Tij{z), = 0, (2-7) 

and from V^^F'"'' = 0, 

= V,^"^ = i-e'^''Z),(r"So), (2-8a) 

= ^^^" = Djr^. (2-8b) 
These equations imply that the electric field Eq takes the Coulomb form, 

^0 = ^, (2-9) 

and = ^J^ij{z)dz'' A dz^ is a harmonic form on JC^. In general, there may exist 
such a harmonic form that produces an energy-momentum tensor consistent with 
the structure of the Ricci tensors in H2-5() . if the second Betti number of /C^ is not 
zero. The monopole-type magnetic field in four-dimensional spacetime provides such 
an example. However, since scalar and vector perturbations become coupled if such 
a background field exists, in the present paper we consider only the case Tij = 0. 

With this assumption, the energy-momentum tensor for the electromagnetic 
field, ^ 

— J' 'J "^jjLv-F '^fjJ' ^ 1 (2'10) 

is written 

= -P6l Ti = P5i- P = \eI = (2-11) 
and the background Einstein equations, 

Rf^v = ^Ag^^ + (t^^ - l^Tg^A , (2-12) 



are reduced to 



where 



2(n — 

nnr = r'^R-2{n + l)Xr + ^ 2n-i ' (2-13a) 
2DaDbr = arg^b, (2-13c) 



A:=^^, Q'--=^^y (2-14) 

n(n -|- 1) n[n — 1) 
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From (|2-13a;) , (|2-13c;) and the identity 



DaUr = UDaV - ^^RDar, (2-15) 



it follows that 

j.n—1 

Hence, we obtain 



0. (2-16) 



, , 2n(n-l)M 2(n - l)(2n - 1)0^ 

^i? = 2A + ^ ^ ^^T^ —. (2-17c) 

When Vr 7^ 0, these equations give the black hole type solution 

ds^ = -f{r)dt' + + r'^dal, (2-18) 

with 

/W=i^-Ar^-^ + ^. (2-19) 

To be precise, the spacetime described by this metric contains a regular black hole for 
Ac < A < if K = or = -1, and for A^i < A < Xc2 and Q'^/M'^ <{n + lf/{^n) 
\i K = 1, where Ac, Ad and Ac2 are functions oi Q,M and n. (For details, see 
Appendix 1X1 1 

Next, let us consider the case in which r = a is constant. Here, we obtain the 
Nariai-type solutioii^'l^ 

dn^ 

ds" = -f{p)dt^ '^JIp)^ ^^'^^^ 



where 

o = \n -^ ijA 

a 

and the constant a is determined solution to 



/(p) = l-ap2; ^ = (n + l)A- Jf^' , (2-21) 



(n-l)^ = (n + l)A+ ^" (2-22) 

2.2. Tensorial decomposition of perturbations and the Einstein equations 

In general, as tensors on /C^, the metric perturbation variables h^j^y = Sg^u are 
classified into three groups of components, the scalar hat, the vector hai and the 
tensor hij. Unfortunately, this grouping is not so useful, since components belonging 
to different groups are coupled through contraction with the metric tensor and the 
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covariant derivatives in the Einstein equations. However, in the case that /C" is 
maximaUy symmetric, if we further decompose the vector and tensor as 



ha^ = Diha + h]^/ ; &h]^/ = 0, (2-23a) 
hij = hLjij + hTij] hr] = 0, (2-23b) 

D^h^^ = 0, D^h% = 0, (2-23c) 



the Einstein equations are decomposed into three groups, each of which contains only 
variables belonging to one of the three sets of variables {hab, ha,hL, h^}, {h^^^ , hipj} 

and {hPj}W^^ Variables belonging to each set are called the scalar-type, the 
vector-type and the tensor-type variables, respectively. This phenomenon arises be- 
cause the metric tensor 'jij is the only non-trivial tensor in the maximally symmetric 
space, and as a consequence, the tensorial operations on h^u to construct the Einstein 
tensors preserve this decomposition.^S Moreover, for the same reason, the covariant 
derivatives are always combined into the Laplacian in the Einstein equations after 
this decomposition. Thus, the harmonic expansion of the perturbation variables with 
respect to the Laplacian is useful. 

In the case in which /C" is of the Einstein type, the Laplacian preserves the 
transverse condition, 

DW ■ bhlp. = 0, (2-24) 

and (|2-23|) leads to the relations 

b'hai = Aha, (2-25a) 
D^hTij = ^^^^bi{A + nK)hP + [b-b + {n- l)K^ h^^, (2-25b) 

D'b^hTij = ^^^^(^ + nK)AhP. (2-25c) 

Hence, the tensorial decomposition (|2-23|) is still well-defined if these equations can 
be solved with respect to ha, h^ and h^^. We assume that this condition is satisfied 
in the present paper. 

The relations ()2-25|) also show that tensor operations on h^^, that lower the 
rank as tensors on /C" preserve the tensorial decomposition into the scalar type and 
vector type, because the Weyl tensor C^jki of IC^ is of second order with respect to 
differentiation and does not take part in such operations. It is also clear that tensor 
operations on hab and hai that preserve or increase the rank have the same property. 
Furthermore, the covariant derivatives are always combined into the Laplacian in 
the tensor equations for the scalar-type and vector-type variables obtained through 
these operations. Therefore the difference between the perturbation equations in 
the Einstein case and the maximally symmetric case can arise only through the 
operations that produce the second-rank terms in {SGij)T from hij. As shown by 
Gibbons and Hartnoll,^ these terms are given by 

ALhij := -b ■ bhij - 2Rikjih''^ + 2(n - l)Khij. (2-26) 
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The Lichnerowicz operator A l defined by this equation preserves the transverse and 
trace- free property of hij as pointed out in Ref. i36)l . Furthermore, we can easily 
check that the following relations hold: 

Kl {o^bj - ^7*jA^ 4°^ = - {bib J - ^7ijA^ KhP, (2-27a) 
Ai ibuh^',],) = -bu (b-b-{n- 1)k) h^'l (2.27b) 



Hence, the Lichnerowicz operator also preserves the tensorial types. One can also see 
that only the Laplacian appears as a differential operator in the perturbed Einstein 
equations for the scalar-type and vector-type components after the tensorial decom- 
position. Thus, we can utilize the expansion in terms of scalar and vector harmonics 
for these types of perturbations, and the structure of the Einstein space affects the 
perturbation equations only through the spectrum of the Laplacian D ■ b. Similarly, 

(2) 

for tensor-type perturbations, if we expand the perturbation variable /i^/^ in the 
eigenfunctions of the Lichnerowicz operator, we obtain the perturbation equation 
for the Einstein case from that for the maximally symmetric case by replacing the 
eigenvalue A;|, for —b ■ b in the latter case by — 2ni^, where is an eigenvalue 
of Ai. 

2.3. Perturbation of electromagnetic fields 

The Maxwell equations consist of two sets of equations for the electromagnetic 
field strength T = {l/2)T^^dxf' A dx" , dT = and V^J^'"' = Ji". If we regard 
{6T)fj,u = 5!F^y as the basic perturbation variable, the perturbation of the first 
equation does not couple to the metric perturbation, and it is simply given by 

d5T = 0. (2-28) 

This is simply the condition that 6 J- is expressed in terms of the perturbation of the 
vector potential, 5 A, as 5J^ = d5A. Next, perturbation of the second equation gives 
two set of equations, 

^D,{r'\5Tr') + bi{5J^r + Eoe'^' (^^f (/^^ " K) " D,h}}j = J\ (2-29a) 
^ Z)Jr"-2((5^),°+Eoe"%b) + Dj{8F).i:> = J,, (2-29b) 



j,n—2 

where h^y = 6g^u, and the external current is treated as a first-order quantity. 

The perturbation variable 6J- transforms under an infinitesimal coordinate trans- 
formation — > + 5x^ as 

SiSJ'^..) = -{&sJT)^, = -WV^T^. + 2.F,[^V,]fe". (2-30) 

To be explicit, writing bx^ as 



(2-31) 
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we have 

5{^:Fab) = -D,{EoT'')eab, (2-32a) 

6{^J^a^) = -EoeabD^T', (2-32b) 

6{5J='ij) = 0. (2-32c) 

As in the case of the metric perturbation, the perturbation of the electromagnetic 
field, 5J-^y, can be decomposed into different tensorial types. The only difference is 
that no tensor perturbation exists for the electromagnetic field, because it can be 
described by a vector potential. After harmonic expansion, the Maxwell equations 
(|2-29j) are decomposed into decoupled gauge-invariant equations for each type in the 
background consisting of H2-18() and H2-2U() . because the Weyl tensor of /C" does not 
appear in the Maxwell equations. This gauge- invariant formulation is given in the 
subsequent sections for each type. Here, we only note that is gauge-invariant for 
a vector perturbation, since from H2-32() the gauge transformation of 5J- does not 
depend on U. In contrast, for a scalar perturbation, these perturbation variables 
must be combined with perturbation variables for the metric to construct a basis for 
gauge-invariant variables for the electromagnetic fields. 

Finally, we give expressions for the contribution of electromagnetic field pertur- 
bations to the energy-momentum tensor: 

5Tab = ^ {e'''6Tcd + Eohij gab - ^Eihab, (2-33a) 
6Tt = E^e'^HT.b, (2-33b) 
m= ie^'Sr^d + Eohl) 6). (2.33c) 



§3. Tensor-type Perturbation 



Because an electromagnetic field perturbation does not have a tensor-type com- 
ponent, the electromagnetic fields enter the equations for a tensor perturbation only 
through their effect on the background geometry. 

As explained in H2.21 tensor perturbations of the metric and the energy- momentum 
tensor can be expanded in terms of the eigentensors Tj-,- of the Lichnerowicz operator 
Al defined by (t^^ as 

Sgab = 0, 5gai = 0, Sgij = 2r^HT{y)Tij, (3-la) 
5Tab = 0, ST,^ = 0, 6T; = TT(y)T;.. (3-lb) 

The expansion coefficients Ht and tt themselves are gauge- invariant, and the Ein- 
stein equations for them are obtained from those in the maximally symmetric case 
considered in Ref. 141(1 through the replacement k'^ — > — 2nK. The result is 
expressed by the single equation 

HHt + ^Dr ■ DHt - ^''"^^^^"^^^ ffr = -k\t, (3-2) 
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where A/, is the eigenvalue of the Lichnerowicz operator, 



(3-3) 



For the Nariai-type background H2-2U|) . the above perturbation equation simph- 
fies to 



-dfHT + fdpifdpHr) 



Xl - 2(n - 1)K 



For the black hole background (|2-18|) . if we introduce the master variable ^ by 
(|3-2() can be put into the canonical form 



T 



where 



Vt 



f 



nr f n(n — 2)f 



2 4 

In particular, for /(r) given by H2-19|) . Vt is expressed as 

- lOn + 8 _ n{n + 2) ^ n{3n - 2)Q^ 



Vt = ^ 



Xl + 



-K 



-Xr' + 



2r 



n-l 



(3-4) 
(3-5) 

(3-6) 
(3-7) 

(3-8) 



Note that 1)3 -Gf) with ()3-7() describes the behaviour of a tensor perturbation if 
the background metric has the form H2-18() . Hence, it may apply to a system more 
general than the Einstein-Maxwell system considered in the present paper. 



§4. Vector-type Perturbation 

We expand vector perturbations in vector harmonics Vj satisfying 
{D-D + kl)Yi = 0; b'Yi = 
and the symmetric trace-free tensors Vjj defined by 



Note that this tensor is an eigentensor of the Lichnerowicz operator, 

AiV,, = [kl + {n-l)K]Nij, 



(4-1) 



(4-2) 



(4-3) 



but it is not an eigentensor of the Laplacian in general. 

In this paper, we assume that the Laplacian —D-D is extended to a non-negative 
self-adjoint operator in the L2-space of divergence-free vector fields on /C", in order 
to guarantee the completeness of the vector harmonics. Because —D-D is symmetric 
and non-negative in the space consisting of smooth divergence-free vector fields with 
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compact support, it always possesses a Friedrichs extension that has the desired 
property.^ With this assumption, ky is non-negative. 

One subtlety that arises in this harmonic expansion concerns the zero modes of 
the Laplacian. If /C" is closed, from the integration of the identity D^[V^ DiVj) — 
D'VWiVj = VW ■ DVj, it follows that DiYj = for k'^ = 0. Hence, we cannot 
construct a harmonic tensor from such a vector harmonic. We obtain the same result 
even in the case in which /C" is open if we require that Y^DiYj fall off sufficiently 
rapidly at infinity. In the present paper, we assume that this fall-off condition is 
satisfied. From the identity D^DiVj = DiD^Vj + {n — l)KVi, such a zero mode 
exists only in the case K = {). We can further show that vector fields satisfying 
DiVj = exist if and only if /C" is a product of a locally flat space and an Einstein 
manifold with vanishing Ricci tensor. 

More generally, Vjj vanishes if Vj is a Killing vector. In this case, from the 
relation 

2kvDjY\ = [kl - (n - l)K] Yi, (4-4) 

ky takes the special value ky = {n — 1)K. Because ky > 0, this occurs only for 
K = Q OT K = \. In the case = 0, this mode corresponds to the zero mode 
discussed above. In the case K = \, since we are assuming that /C" is complete, 
/C" is compact and closed, as known from Myers' theorem,^ and we can show 
the converse, i.e. that if ky = [n — 1)K, then Yij vanishes, by integrating the 
identity = Y'*&Yij = D^{Y'*Yij) + kvY'^*Yij over /C". Furthermore, using the 
same identity for Y' D^Yij, we can show that there is no eigenvalue in the range 
0<kl <n- l.E^ 

4.1. The Einstein equations with a source 

In terms of vector harmonics, perturbations of the metric and the energy- 
momentum tensor can be expanded as 

Sgab = 0, 6gai = rfaYi, 6gij = Ir'^HxYij, (4-5a) 
STab = 0, 6Tt = rr"Vi, ST; = ttY). (4-5b) 

For my = k'^ — {n — 1)K / 0, the matter variables Tq and tt are themselves 
gauge-invariant, and the combination 

i.„ = /„ + .«„(|£) ,4.6) 

can be adopted as a basis for gauge-invariant variables for the metric perturbation. 
The perturbation of the Einstein equations reduces to 

Da (r"+^F(i)) - mvr^-'eabF' = -2/^V"+l6„,r^ (4.7a) 
kvDa{r''-^F'') = -K^r'^TT, (4-7b) 
where Eah is the two-dimensional Levi-Civita tensor for gab, and 

^(1) = ^ab^jj^ ^ ^ab^jj^ (^f^\^ ^ (4.g) 
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This should be supplemented by the perturbation of the energy-momentum conser- 
vation law 

^^(^n+i^a^ + ^r"rT = 0. (4-9) 
2kv 

Now, note that, for my = 0, the perturbation variables Ht and tt do not exist. 
The matter variable is still gauge-invariant, but concerning the metric variables, 
only the combination F^^^ defined in terms of fa in (|4-8|) is gauge invariant. In 
this case, the Einstein equations are reduced to the single equation 1)4 -Zaj) . and the 
energy-momentum conservation law is given by (|4-9|) without the tt term. 

Now, we show that these gauge-invariant perturbation equations can be reduced 
to a single wave equation with a source in the two-dimensional spacetime AA^. We 
first treat the case my ^ 0. In this case, by eliminating tt in H4-7b|) with the help 
of (jHH), we obtain 

Dair'^'^F") = Dair^'+W"). (4-10) 

my 

From this it follows that F"' can be written in terms of a variable as 

^n-lpa ^ ^abj^^^ ^ r'^+V'^. (4-11) 

my 

Inserting this expression into H4-7a() . we obtain the master equation 



r^Da ( ^D'^f^) - = -^r"e"^Z),(rTb). (4-12) 



r" y my 



Next, we consider the special modes with my = 0. For these modes, from (|4-9|1 
with Tt = 0, it follows that can be expressed in terms of a function r(i) as 

r"+Va = e,,DV«. (4-13) 

Inserting this expression into (|^^ with e'"^Dc{Fd/r) replaced by F^'^'^/r, we obtain 

D„(r"+^FW) = -2K2L)„rW. (4-14) 

Taking account of the freedom of adding a constant in the definition of t^^\ the 
general solution can be written 

Hence, there exists no dynamical freedom in these special modes. In particular, in 
the source- free case in which r^^^ is a constant and K = 1, this solution corresponds 
to adding a small rotation to the background solution. 

4.2. Einstein- Maxwell system 

As shown in ^2.3( all components of SJ^^u are invariant under a coordinate gauge 
transformation for a vector perturbation. In order to find an independent gauge- 
invariant variable, we expand 5J^ai in vector harmonics as 6J^ai = ^aVj. Then, since 
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^^ab = for a vector perturbation, the (a, h, i)-component of the Maxwell equation 
((T^ is written 

D^aA] = 0. (4-16) 

Hence, Aa can be expressed in terms of a function A as Aa = DaA. Then, the 
(a, i, j)-component of (|2-28() is written 

DaSJ'ij = DaA{d^Vj - djYi). (4-17) 

This implies that 6J^ij can be expressed as 

5J^,j = + A{d,Yj - djY,), (4-18) 

where Cij is an antisymmetric tensor on /C" that does not depend on the y-coordinates. 
Finally, from the A;)-component of 1)2-28^ . it follows that Cijdz'^ A dz^ is a closed 
form on /C". Hence, Cij can be expressed in terms of a divergence- free vector field 
Wi as Cij = diWj — djWi- With the expansion in vector harmonics, we can assume 
that Wi is a constant multiple of Vj, without loss of generality. Hence, this term can 
be absorbed into A by redefining it through the addition of a constant. Thus, we 
find a vector perturbation of the electromagnetic field can be expressed in terms of 
the single gauge-invariant variable A as 

5Fab = 0, 5Fai = DaAY,, bT,, = A (AVj - DjYi) . (4-19) 

Next, we express the remaining Maxwell equations in terms of this gauge- 
invariant variable. For a vector perturbation, only Ea. H2-29b|) is non-trivial. If we 
expand the current Jj as 

Ji = JV„ (4-20) 

it can be written 

(r"-2i?"^ - r"-ii?oe"Vfe) + AD^ (aV,- - DjYi^ = JYi. (4-21) 
Hence, from the identity 

(Av, - A-v.) = 4 + iy m ^^^ (4.22) 

the gauge-invariant form for the Maxwell equation is given by 

-L_Da{r^^-'D-A) - Ml±i!izl)^^ = _ J + rEoF^'l (4-23) 

In order to complete the formulation of the basic perturbation equations, we 
must separate the contribution of the electromagnetic field to the source term in the 
Einstein equation (|4-12() . Because for a vector perturbation, (|2-3c{|) is expressed as 

6T^ = 0, 6T^^ = -Eoe^'^Di^AY^, ST] = 0, (4-24) 

the contributions of the electromagnetic field to and tt are given by 

= -^e,,Z)^^, 4^ = ^. (4.25) 
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Hence, the Einstein equations for the Einstein-Maxwell system can be obtained by 
replacing Tq in (|4-12|1 by 

ra = r^^ + fa, (4-26) 

where the second term represents the contribution from matter other than the elec- 
tromagnetic field. 

In order to rewrite the basic equations obtained to this point in simpler forms, 
we treat the generic modes and the exceptional modes separately. 

4.2.1. Generic modes 

For modes with mv = ky — {n — 1)K ^ 0, the above simple replacement yields 
a wave equation for with oA in the source term. This second derivative term can 
be eliminated if we extract the contribution of the electromagnetic field to as 

n = n- (4-27) 

my 

The result is 

r^Da (^D-f2) -^f2 = ^A - ^r\-^Da{rn). (4-28) 
Equations (|4-11() and H4-9p are replaced by 

^n-l^a ^ ^abjj^^ ^ r"+^f ^ (4-29a) 

my 

Dair^'+^n + ^r^TT = 0. (4-29b) 
2kv 

Finally, inserting this expression for into (|4-23|) and using (|4-28j) . we obtain 
^.D^ir^-^D'^A) -^j4 + in- l)K + ^^^^^^) A='^n-J. (4-30) 



Thus, the coupled system of equations consisting of (|4-28|) and (|4-3()|) provides the 

"V 



basic gauge-invariant equations with source for a vector perturbation with ky ^ 



{n-l)K. 

4.2.2. Exceptional modes 

For exceptional modes with ky = {n — 1)K, from the definition H4-13() of t^^^ 
and (|4-25() . t^^^ can be expressed as 

rW = -qA + f^^\ (4-31) 

Hence, ()4-15() can be rewritten as 

Inserting this into (|4-23j) . we obtain 

_i,Z>„(...-Z,M) - J, (2(„ - + ^ = _ |^,<.). ,4.33) 

Therefore, only the electromagnetic perturbation is dynamical. 
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4.3. Decoupled master equations 

The basic equations for the generic modes are coupled differential equations and 
are not useful. Fortunately, they can be transformed into a set of two decoupled 
equations by simply introducing master variables written as linear combinations of 
the original gauge-invariant variables. Such combinations can be found through 
simple algebraic manipulations. 

4.3.1. Black hole background 

For the black hole background (|2-18|) . appropriate combinations are given by 

^± = a±r-"/2/2 + 6±r"/2-i^, (4.34) 

with 

^"+'^+^=( (n^-'^lTM + Z^ 'f)' 

where Z\ is a positive constant satisfying 

A^ = {n'^-lfM'^ + 2n{n-l)mvQ^. (4-36) 

When expressed in terms of these master variables, H4-28() and H4-3U|) are transformed 
into the two decoupled wave equations 

- -7^^± = Sv±. (4-37) 



Here, 



Vv± = ^ 



where 



2 (n^ -2n + 4)K n{n - 2) 2 ^(5"- - '^)Q A^i 
1 1 4^2n-2 ^ 



(4-38) 



and 



li± = ± Z\, (4-39) 

5y± = -a± -e'^^'Dairn) - b^r^'^-^fj. (4-40) 

my 

For n = 2,K = I and A = 0, the variables and are proportional to the 
variables for the axial modes, z[ ^ and Z2 ^ given in Ref. IHS)) . and Vy+ and Vy- 

coincide with the corresponding potentials, ^ and V2 \ respectively. 

Here, note that in the limit Q — > 0, ^+ becomes proportional to A and to 
Hence, and represent the electromagnetic mode and the gravitational mode, 
respectively. In particular, in the limit Q ^ 0, the equation for ^_ coincides with 
the master equation for a vector perturbation on a neutral black hole background 
derived in Paper I. 
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4.3.2. Nariai-type background 

For the Nariai-type background (|2-2U() . the combinations 

<P-i- = a±fi + b±A, 

with 

mvQ Q 



(a_,6_)= 1, 



Q2n-2 [2{n - 1)K - o?a + a? An] ' q 
2n{n - l)g2 



q[2{n - 1)K - a^a + u'^An] 



(4-41) 

(4-42a) 
(4-42b) 



give the decoupled equations 



4 + (n - 1)K 



a ± An <P± 



2^ 

my 



a±fa^+^e'>^D,f, - b±fJ, 



where 



A 



N 



2 , 2n(n-l)Q2 



C7^ + 



{ey + {n-l)K} 



1/2 



(4-43) 
(4-44) 



§5. Scalar-type Perturbation 

We expand scalar perturbations in scalar harmonics satisfying 



(5-1) 



As in the case of vector harmonics, we assume that —A is extended to a non- 
negative self-adjoint operator in the L^-space of functions on /C". Hence, k"^ > 
0. In the present case, such an extension is unique, since we are assuming that 
/C" is complete,!^ and it is given by the Friedrichs self-adjoint extension of the 
symmetric and non-negative operator — A on C^{K7^). If /C" is closed, the spectrum 
of A is completely discrete, each eigenvalue has a finite multiplicity, and the lowest 
eigenvalue is /c^=0, whose eigenfunction is a constant. A perturbation corresponding 
to such a constant mode represents a variation of the parameters A, M, and Q of the 
background, as seen from the argument in ^2.11 For this reason, we do not consider 
the modes with k"^ = 0, although there may be a non-trivial eigenfunction with 
A;^ = in the cases in which /C" is open. Note that when A;^ = is contained in the 
full spectrum but does not belong to the point spectrum, as in the case IC^ = M", it 
can be ignored without loss of generality. 

For modes with A;^ > 0, we can use the vector fields and the symmetric trace-free 
tensor fields defined by 



(5-2a) 
(5-2b) 
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to expand vector and symmetric trace-free tensor fields, respectively. Note that Sj 
is also an eigenmode of the operator D ■ D, i.e., 

[b-b + k^ -{n- l)K]^^ = 0, (5-3) 

while Sjj is an eigenmode of the Lichnerowicz operator: 

(Ai - k')^,j = 0. (5-4) 

In the case of scalar harmonics, the modes with k"^ = nK are exceptional. Given 
our assumption, these modes exist only for = 1. Because /C" is compact and 
closed in this case, from the identity 

m = ^"S.. (5.5) 

we have bj^l = 0. From this, it follows that / dJ^ z^/^§ij*S'^^ = 0. Hence, §ij 
vanishes identically. 

5.1. Metric perturbations 

In terms of scalar harmonics, perturbations of the metric and the energy-momentum 
tensor are expanded as 

Sgab = /afeS, Sgai = rfa^i, 6gij = 2r^{HLjijS + HTSij), (5-6a) 
6Tab = TabS, 5T^ = rr'^Si, 5T} = 5P6p + tt^. (5-6b) 

Following Ref. 141(1 . we adopt the following combinations of these expansion coeffi- 
cients as the basic gauge-invariant variables for perturbations of the metric and the 
energy-momentum tensor: 

F = Hl + -Ht + -D^rXa, (5-7a) 
n r 

Fab = fab + DaXb + DbXa, (5-7b) 

^ab = Tab - P{DaXb + DbXa) - X'D.Pgab, (5-7c) 

2k 

Sa = Ta + —PXa, (5-7d) 



r 



Ul = 6P + X^'DaP. (5-7e) 



Here, 



= ^ [fa + ^DaHr) , (5-8) 

and we have used the background value for T^^ given in H2-lip . 

Now, note that for the exceptional modes with k"^ = n for K = 1, Ht and tt are 
not defined, because a second-rank symmetric tensor cannot be constructed from §, 
as mentioned above. For such a mode, we define F, Fab, ^abi and Ul by setting 
Ht = in the above definitions. These quantities defined in this way are, however, 
gauge dependent. These exceptional modes are treated in Appendix IdI 
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5.2. Maxwell equations 

Because Xa defined in (|5-8|) transforms under (|2-,S1|) as 

SXa = Ta, (5-9) 

from ()2-32|) it follows that the following combinations £ and £a can be used as basic 
gauge-invariants for perturbations of the electromagnetic field: 

6J^ab + D,{EoX'')€abS = £eabS, (5-lOa) 
6Tai - kEoeabX% = reab£%. (5- 10b) 

Here, note that 6J^ij = for a scalar perturbation, since 6J-ij can be written as 
6Tij = Di6Aj — Dj6Ai, from the Maxwell equations, and SAi is just the gradient of 
a scalar perturbation. 

By expanding the current Jj as 

Ji = rJSi, (5-11) 
the Maxwell equations ()2-29|) can be written 

^Dair^S) + -£a - ^Da{F^ - 2nF) = eabJ\ (5-12a) 
e''^Da{r''-^£b) = -r''~^J, (5- 12b) 



while (|2-28|) is expressed as 

£ = -'^D,{r£^). (5-13) 

Note that (|5-12a|) and (|5-12bj) give the current conservation law 

Dcir'^r) = -kr^'-^J. (5-14) 

We can reduce these equations to a single wave equation for a single master 
variable. First, from ()5-12b|) and 1)5- 14() . we have 

e''''Da{r^-'£b) = jD.ir'^r). (5-15) 
k 

Therefore, if we define Ja by 

r = ^e^'Jb, (5-16) 
£a can be expressed in terms of a function A as 

£a = :^[DaA + Ja). (5-17) 

Then, the insertion of this into (|5-12a() yields 

k^D.A = -Da{r^£) + lDa{F^^ - 2nF). (5-18) 
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Therefore (adding some constant to A in its definition, if necessary), we obtain 

= -k'^A + ^{F^ - 2nF). (5-19) 

Thus, the gauge-invariant variables £■ and £-a can be expressed in terms of the single 
master variable A. Finally, by inserting these expressions into (|5-13p . we obtain the 
following wave equation for A: 

r-'D^ (^) - = -r-'l^a (^^) - ^(F.^ - 2nF). (5-20) 

Next, we derive an expression in terms of A for the contribution of the electro- 
magnetic field to the perturbation of the energy- momentum tensor. First, for a scalar 
perturbation, H2-33() can be written in terms of £, Ea and the metric perturbation 
variables as 

r,fe = E^\~£^- D.iEoX'^) + ^Eof^^ gab - \Elfai>. (5-21a) 

+ ^E^X'^^ , (5-21b) 

D,{EoX') - ii?o/c) , (5-21c) 

TTT = 0. (5-21d) 

Inserting these into (|5-7c|) - (|5-7e|) and using (|5-17|) and (|5-19j) . we obtain the following 
expressions for the corresponding gauge-invariant variables Eat, and El in terms 
of A, F and Fab'- 



Tab = 


Eo 


t" = 


-^0 ( 


SP = 


Eo ( ' 



=[^-^A + ^F]gab- ^Fab, (5-22a) 



S!lA+^F]q,-^ 
^r = -^,{DaA + Ja), (5.22b) 

5.3. Einstein-Maxwell system: Black hole background 

For generic modes of scalar perturbations, the Einstein equations consist of four 
sets of equations of the forms 

Eab = K^Eab, Ea = Ea, El = K^Ul, Et = K^T- (5-23) 

(For the definitions of Eab,Ea,EL and Ft, see Eqs. (63)-(66) in RefHTJ.) If we 
introduce the variables F^^ and F as 

the equation for Ft is algebraic and can be written 

2{n - 2)F + F^ = -St, (5-25) 
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where 

St = -p-'^^^T. (5-26) 
Hence, if we introduce X, Y and Z defined by 

X = Fl- 2F, Y = F^ - 2F, Z = F[, (5-27) 
as in Paper I, the original variables are expressed as 

F = -^{X + Y + St), (5-28a) 
2n 

Fl = X + 2F, F; = Y + 2F, F[ = Z. (5-28b) 

In contrast, for the exceptional modes, (|5-25)) is not obtained from the Einstein 
equations. However, this equation with St = can be imposed as a gauge condition, 
as shown in Paper I. Under this gauge condition, all equations derived in this sec- 
tion hold without change. However, the variables still contain some residual gauge 
freedom, and we must eliminate it in order to extract physical degrees of freedom. 
This is done in Appendix IdI 

To obtain the expressions of the remaining Einstein equations in terms of X, Y 
and Z, we introduce E^, Ea and E^ defined by 

Eab = r'^-^Eab + ^„+2 {Fab " 2nFgab), (5-29a) 
Ea = r'^-^Ea, (5-29b) 
E, = r^-^E, - ^^^^^^F, (5.29c) 

and separate the contributions of the electromagnetic field to the perturbation of the 
energy-momentum tensor as 

n-l 2 

Sab — r l^[^ab-^ab ),^a — [^a - ^a J^Sl — V K[2jL-^L )■ 

(5-30) 

Then, the remaining Einstein equations are written 

^"■b = ~Sab + Sab, (5-31a) 

. kn{n-l)Q^ DaA + Ja , 

^« = ^^in — + (^-^ib) 

^L = ^2 — 7 + ^^- (^'^1^) 

As shown in Paper I, if the equations corresponding to Ea, El and El are 
satisfied, the other equations are automatically satisfied, as seen from the Bianchi 
identities, provided that the energy momentum tensor satisfies the conservation law. 
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which in the present case is given by the two equations 



^Da{r S )-Sl + ^ St - 0, 

^Db{r SJ + ^Sa - u-^Sl = k —^Ja- 



(5-32a) 
(5-32b) 



The exphcit expressions of the relevant equations in terms of X, Y and Z are as 
fohows. First, the parts and the part are 



{2r/k)Et = -dtX - drZ - dtSr 
{2r/k)Er ■- 



f 1 77 — 1 

drY +J-(X-Y) + -^dtZ - drSr + St, 

2/ r 



(5-33a) 
(5-33b) 



2EI = fdtdrX - ( ^/ + t]a^X + fdtdrY + 



2/ /' 



dfY 



+ ^Z + fdtdrST 



/ + y ) dtST. 



(5-33c) 



After the Fourier transformation with respect to the time coordinate t, setting 
dt{X,Y,Z) = -iuj{X,Y,Z) and solving with respect to {X',Y',Z') = dr{X,Y,Z), 
we obtain 



, n- 2 
X' = X + 



f 2 



Y + 



-2k'^Eo (^M + J^^ + ( l-^ --^ St + 2Sr 



2/ r 



2/ luj 



+2k^Eq [a! + Jr^ -S't + "^—^St - 



2Sr, 



Z' 



ILO 



X - 2k'^Eq a 



Jt 



+ St 



2St 



(5-34a) 

(5-34b) 

(5-34c) 



Next, El is expressed as 



1 



2EI = -dfx - '—drX + -dfY -\ — + —\drY + —dtZ 



f 

+ 
+ 



/ 

(n-l)(n-2)M 



/' , nf 



2n. 



- (2n- 1)A + 



2 ' r J ' rf 
(n-l)Q2 



r.2n 



k'^-nK {n+l)(n-2)M 



1 

7^ 



+-,dfST -['- + 



r-n+1 



drST 



+ (n + 1)A + 



X 

(2n^ - 2n - 1)Q^ 



Y 



^k^ + n{n-2)K (n + l)(n-2)M 



+ 



r-n+l 



21 



(n2 - 1)A + 



(n^ -n + 1)Q 



21 



^2n 



St. 



(5-35) 



Applying the Fourier transformation with respect to time to the corresponding equa- 
tion and eliminating X' , Y' and Z', with the help of (|5-34j) . we obtain the following 
linear constraint on X, Y and Z: 

2 n{n-l)M (n - l)(2n - 1)Q2 



M 



„n—l 



-.2 



2n-2 



—^[n{n + l)\r' + 



(n^-l)M 3n(n-l)Q^ 



r,n—l 



+ 



r2n-2 



LO^r'^ - k^f + nK'^ - (n - l)KXr 



n{2n — l)Ar 



2 n(n — 1)(5^ 



r,2n-2 



2 AKnM K{n^ - An + 1)Q2 



r.n— 1 



M ({n + lfM AnQ"^ 



n—l 



r2n-2 



+ 



,2n-2 



n{n — l)Ar^ + 



^2n-2 



r.2n-2 



{Y + St) 



nuj'^r'^ + ( Xr 



2 (n-l)M , (n-l)Q2^^^2 



„n— 1 



+ 



„2n-2 



nrf{A' + Jr) + k^A 



-2nrfSr - r^f— + 2r^fS: = 0. 
iu> 



(5-36) 



Using almost the same method as in Paper I, we can reduce this constrained 
system for X, Y and Z to a single second-order ODE with source. The master 
variable in the present case is given by 



n 



■i/2-2fJ 



rpr 

-^ + 2F 

iujr 



X + Y + St-^ 

j,n/2-2jJ 



where 



n(n+l) 2 
H = m -\ ^ X — n z, 

m = k^ — nK, 



2M 



X 



^, y = Xr , z 



r.2n-2 ■ 



(5-37) 

(5-38) 
(5-39) 
(5-40) 



This master variable coincides with that for the neutral and source-free case in Paper 
I for Q = and St = 0. The master equation for (p is 



fif^y + {uj^ - vs)^ = s^. 

Here, the effective potential Vs is expressed as 
f{r)Us{r) 



Vs 



16^2^2 ' 



(5-41) 



(5-42) 



22 



Us = [-n^{n + 2)(n + l)^x^ + 4n^(n + l){n{r? + 6n - i)z + 3(n - 2)m}x 
-12n^(3n - 2)2^ - 871^(1171^ - 26n + 12)mz - 4(n - 2)(n - 4)™^] y 
+n'^{n + + n(n + 1) {-3n^(5n^ - 5n + 2)z + 4(2n^ - 3n + 4)m 

- 2)(n - 4)(n + 1)K} + 4n [n^(4n^ + 5n^ - lOn + 4)z^ 
- {n(34 - 43n + 21n^)m + n^in + l)(n^ - lOn + 12)K} z 
-3(n - A)m^ - 3n(n + l)(n - 2)Km\ x - 4n^(3n - 2)z^ 



+I2n^ {2(-6n + 371^ + 4)m + n^(3n - 4)(n - 2)K} z 
+ {4(13n - 4)(n - 2)m^ + Sn^illn^ - 18n + A)Km} 
+16m^ + 4n(n + 2)Km^, 

and the source term S,p has the fohowing structure: 



(5-43) 



Sip 



f 



Psi 
H 



A 



Jt 



2 Jt 



+ 2nrfJr + 2k^— + 2nf 



,rdrJt 



UjJ 



.HSt - - 2n/^ - 2nrfSr 



H IUJ 



H iu) 



(5-44) 



where Psi,Ps2 and P53 are functions of r expressed as polynomials of x,y and z. 
Their explicit expressions are given in Appendix O 

The basic variables X, Y and Z are expressed in terms of the master variable ^ 

as 



X = r"/2-2 



lujr 



n/2-1 





Pxo \ 


, / 


mm ) 




Pyo 


/ 






frdr<P 



Px\ 
AH 

+ Zs, 



(5-45a) 
(5-45b) 
(5-45c) 



where XgyYg and Zg are contributions from the source terms given by 



2/^2' 



1; 



.riPx2S^,2rvrl , nPx^rSl 2r\^ 
^ 2m iuj^ H m ioj H^^' 

2 „ / Pya . , 2nr/ . „ , , j . nPys ^ 

nPy2 St 2nr/ ^ ^ TiPya rS'[ ^ 2r2 



(5-46a) 



+ TrT77"T— 77— Sr H TTT"^^ H ~ St, 



2m 

2 ^ 2nr/ ^ . ^ , T ^ 2nr/ ^ 2r 



i7 ' ' m iuj 
-iojA + Jf) 



(5-46b) 
(5-46c) 
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Here, Pxo, Pxi, Pyo, Pyi, Pz , PxA, Px2, Pya and Py2 are functions of r and are 
given in Appendix O 

In particular, F is expressed as 



+ —dr^ 

8nH 2 



Hence, using the relation 

- 2nF = -4(n - 1)F - St, 
we can rewrite the Maxwell equation H5-2U() as 



2iujr 



(5-47) 



(5-48) 



2n^{n-l)^zf 



H 



A 



{n-l)q UH^-nPz 
r"/2+2 i 4^ 



<P + frdr^ + 



Sa 



(5-49) 



where 



Sa 



2n\n-lfzp 

7W 



fJr \^2{n-l)Eo 



x.n—2 



iujH 



f{nfSt-rSl). 



(5-50) 

Thus, the two coupled second-order ODEs (|5-41|) and (|5-49|) represent the master 
equations with a source for a scalar perturbation of the Einstein-Maxwell system in 
the background 

Finally, we note that in terms of the variable 



(5-51) 



121) can be rewritten as 
1 



{2H -nrf')Q + 2nrDr ■ DQ - —-{-nfSt + rS^), (5-52a) 



iujH 



1 



where 
and 



Fab + {n- 2)Fgab = - ( DaDbn - -□f?5a6 ]+S^, 



Sa = Sa-K^Eo{DaA+Ja) 



Q 2nr / ~ ~ 2DarDbr ~ 
Sab = {F>arbb + DbVSa J Verb 

-—nfrH (gab - -.DarDbA e'^'^D^ir^ Sa) 



(5-52b) 
(5-53) 



2r r 

+ ^Sab - 



Sc + -St 1 9ab- 



(5-54) 
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These equations are formal, since they contain the factor i/w, which is equivalent 
to integration with respect to the time coordinate t. However, when Sa = K^E^Ja 
and SI = 0, this factor can be eliminated through the replacements St — > wk^EqA 



and e^^Dair^Sb 



-iuj[n 



2)K^r^£'ow4. In this case, H5-52|) provides expressions 



that are manifestly covariant as equations in the 2-dimensional spacetime AA^. The 
source term S$ for <P can also be rewritten in such a covariant form with the help of 
the energy-momentum conservation laws (|5-|-{2|) : 



/ 



r2„n~l ab i 



i 2{nH + 2rH') 



Sb 



+ 



UJ 



rSl + 2r^S^ - HSt 



(5-55) 



The comment made above regarding covariance applies to this expression as well as 

to ^EM- 

5.4. Reduction to decoupled equations 

As in the vector perturbation case, we can transform the master equations (|5-41() 
and (|5-49j) into two decoupled second-order ODEs by introducing new master vari- 
ables written as linear combinations of ^ and A. The only difference in this case is 
that the coefficients are not constant. These new variables are given by 



^± = a±'P + b±A, 



with 



(a+,6^ 
(a-,b- 



(n + l)(M + ^) 
where is a positive constant satisfying 



m {n + l){M + n) (n+l)(M + //)Q 
2r"-i ^' qrn/2-i 
2nQ2 4rig2 



qr 



n/2~l 



, o 4mQ^ 



(n + iy 

The reduced master equations have the canonical forms 
fif^'J + {lo^ - Vs±) ^± = Ss±. 
Here, if we define the parameter 5 by 

H = {1 + 2m6)M, 
the effective potentials Vs± are expressed as 



Vs± 



(5-56) 

(5-57a) 
(5-57b) 

(5-58) 

(5-59) 
(5-60) 

(5-61) 
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with 



H+ = l- !!^!l±}lsx, (5-62a) 
H- = m+ '^'''^^^\ l + m6)x, (5-62b) 



and 



U+ = [-4n^{n + 2){n + ifd'^x'^ - 48n^(n + l)(n - 2)5x 
— 16(n — 2)(n — 4)]y — 6^n'^ 



+45^n^(n + 1)^ {{n + l)(3n - 2)m6 + 4n^ + n - 2} x 



2 



+45(n + 1) {(n - 2)(n - 4)(n + l)(m + n^Er)^ - Tn^ + Tn^ - 14n + 8} 
+ {I6(n + 1) (-4m + 3n^{n - 2)K) 5 - 16(3n - 2)(n - 2)} x 
+64m + 16n(n + 2)K, (5-63a) 
^7- = [-4n^(n + 2)(n + 1)^(1 + mdfx^ + 48n^(n + l)(n - 2)m(l + m5)x 
-16(n - 2)(n - 4)^2] y - n^{3n - 2){n + 1)'^6{1 + m6fx^ 
-4r?{n + 1)^(1 + rabf{(n + l)(3n - 2)mb - r?] x^ 
+4(n + 1)(1 + m6) {m{n - 2)(n - 4)(n + l)(m + n^K)5 
+4n(2n^ - 3n + 4)m + 71^(11 - 2)(n - 4)(n + 1)K} x^ 
-16m { (n + l)m (-4m + 3n^(n - 2)Er) 6 
+3n(n - 4)m + 3n^(n + l)(n - 2)if} x 

+64m^ + 16n(n + 2)m^K. (5-63b) 

The source terms S'si are hnear combinations of S,p = 5*^1^=0 given in H5-44p and 
Sa given in (|5-5()|) : 

Ss± = a±S^ + b±SA- (5-64) 
Here, note that the following relations hold: 

= {n + lfM^6{l + m5), (5-65) 
H = H+H^. (5-66) 

From these relations, we find that Q = corresponds to 5 = 0, and in this limit, 
coincides with ^, and its equation coincides with the mater equation for the master 
variable <P derived in Paper I. Hence, ^_ and represent the gravitational mode 
and the electromagnetic mode, respectively. For n = 2,K = 1 and A = 0, these 
variables and are proportional to the variables for the polar modes, z['^^ 
and ^2^\ appearing in Ref. l35|l . and V5+ and Vs- coincide with the corresponding 
potentials, vj^^ and V2~^\ respectively. 

5.5. Einstein-Maxwell system: Nariai background 

We can derive the master equations for the Einstein-Maxwell system in the 
Nariai-type background H2-2U() in almost the same way as in the black hole case. 
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Actually, the calculations are much simpler in this case. For this reason, we give 
only key equations. As in the previous case, the elimination of the residual gauge 
freedom for the exceptional modes is discussed in Appendix IdI 

First, by defining the variables X, Y and Z by H5-27p . after the Fourier transfor- 
mation, the equations corresponding to Ep, Et and E^ can be written as 



f 

X' = ^Y + 



fa? p J iuj 
f 



2S' 



2k'Eo {A' + Jp]+ i-Sr - T— + 25p, 



2/ 



2/ luj 



+2k^Eo (a' + Jp) -S't- 2S, 



Z' 



^ = X - 2t^EA A - ^ \ + St - 2^ 



UjJ 



lUJ 



The constraint equation obtained from the Ep equation is 



(w^ + a)X + [oj^ + a 



a 



a? ioj 



0. 



These equations can be reduced to a single second-order ODE for 



where 



VF = f[^-2a\, 



Sp 



2na" 



Ji 



2k' Eo 2^ - — -St + 2— + 



St 2a^ a„5f 2a2 



lUJ 



P'^t I op 



UjJ 



Hence, in the present case, we can use F, which is related to X and Y by 

^ ^ X + Y + St 
'~ 2na"-2 ' 

as the master variable. X, Y and Z are expressed in terms of F as 



X 
Y 
Z 



2na" 
2na" 



UJ 



+ a] F- apdpF 



+ a F + opdp 



2^2 

2^^E^A-^SPp, 



fc2 P 



-tJ^ {apF + fdpF) + . 



k^ 



(5-67a) 

(5-67b) 
(5-67c) 



(5-68) 
(5-69) 

(5-70) 
(5-71) 

(5-72) 

(5-73a) 
(5-73b) 
(5-73c) 
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Finally, the Maxwell equation (|5-20|) also holds in the present case, in which it 
becomes 



dp \ dp 



+ uj 



f\A = 2{n- l)a^-'fEoF + 



Sa = IEqSt - iioJt - PdpJp 



(5-74) 



where 

^%Jp- (5-75) 

By taking linear combinations of this equation and H5-69p and introducing the vari- 
ables 

^± = a±F + b±A, 



we obtain the decoupled equations 



dp i/ dp ' ~'~ 



UJ 



cPi = S±, 



with 



S± = a±SF + b^S^A, 
where S,p = S4,\a=o, p is a positive constant satisfying 

4(n-l)2Q2^2 



and 



(a+,6+) 
(a_,6_) 



-,2n+2 



2(n-l)Q Q 

2 + 



2(n- 



p + a, 



qa 



2n 



(5-76) 

(5-77) 
(5-78) 

(5-79) 

(5-80a) 
(5-80b) 



§6. Stability of Black Holes 



In this section, we consider what we are able to conclude at this time about the 
stability of generalised static black holes with charge from the results of our formula- 
tion. In the present paper, we consider only the stability in the static region outside 
the black hole horizon with respect to a perturbation whose support is compact on 
the initial surface. This region is represented as r > for A < and < r < 
for A > 0. As mentioned in ^2.11 such a region exists only for restricted ranges of the 
parameters M, Q and A. These parameter ranges are shown in Fig. ^ (see Appendix 
El for details). 

To study the stability in this region, as in Paper II, we utilize the fact that 
for any perturbation type, the perturbation equations in the static region of the 
spacetime are reduced to an eigenvalue problem of the type 



uj'^<p = A<P, 



(6-1) 
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Fig. 1. The parameter ranges in which the spacetime contains a regular black hole (the shaded 
regions) . 



where ^ is a self-adjoint operator, 

^ = -|| + ^(0; dT, = j, (6-2) 

with V{r) being equal to Vt(?'), yv±{f) or Vs±{r). We regard the black hole to be 
stable if the spectrum of A, i.e., uP', is non- negative. 

To be precise, we must specify a boundary condition for ^ at r = oo in the case 
A < 0, for which the range of r* has an upper bound. In the present paper, we adopt 
the simplest condition, ^ as r ^ oo, in this case, which corresponds to the 
Friedrichs extension of A. For A > 0, the range of is (— oo, +oo) and the operator 
A is essentially self-adjoint; i.e., it has a unique self-adjoint extension.!^ 

In general, if <?(r) is a function with compact support contained m. r > rn (or 
rjy < r < Tc for A > 0), we can rewrite the expectation value of A, (<?, A^)^ as 



dr» 



2 



+ F|^r . (6-3) 



For the Friedrichs extension of A, the lower bound of this quantity coincides with the 
lower bound of the spectrum of A with domain CQ°{r^:). Hence, if we can show that 
the right-hand side of (|6-3() is non-negative, then we can conclude that the system 
is stable. In particular, if V is non-negative, this condition is trivially satisfied. 
However, such a lucky situation is not realized in most cases. One powerful method 
that can be used to show the positivity of A beyond such a simple situation is to 
deform the right-hand side of (|6-3|) by partial integration in terms of a function S as 

{<P, A<P) = J dr^ (^\D<P\'^ + V\<P'^\^ , (64) 
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to (|3-7j) . we obtain 



where 

^ = A + 5, (6-5) 

V = V + f^-S^ (6-6) 
dr 

We call this procedure the S- deformation of V in this paper. As in Paper II, this is 
the main tool for the analysis in the present paper. 

6.1. Tensor perturbation 

If we apply the S'-deformation with 

S = (6-7) 

Vt = ^ [Al - 2(n - l)K] , (6-8) 

irrespective of the r-dependence of /(r). Hence, the effective potential Vp is positive, 
and the system is perturbatively stable with respect to a tensor perturbation if 

Al > 2{n - l)K. (6-9) 

In particular, this guarantees the stability of maximally symmetric black holes for 
K = 1 and K = 0, since A^ is related to the eigenvalue k'^ of the positive operator 
—D • L) as A_L = fcy + 2nK when /C" is maximally symmetric. In contrast, even in 
the maximally symmetric case, Vr becomes negative in the range < /cf, < 2 for 
K = —1, and we cannot conclude anything about the stability for K = —1 from this 
argument alone. 

Note that the condition ()6-9() is just a sufficient condition for stability, and it 
is not a necessary condition in general. In fact, for a tensor perturbation, we can 
obtain stronger stability conditions directly from the positivity of Vr if we restrict 
the range of parameters. For example, for K = 1 and A = 0, it is easy to see that 
Vt is positive if 

AL + 2-2n + ^ ^ > (6-10 

M + v/M2 - g2 - 

2 //'Q^ o^2 



for M^>Q^> 8n{n - l)M7(3n - 2)^, and 

n2 - lOn + 8 , , 

Al + ^ >0 (6-11) 

for < 8n(n — l)M2/(3n — 2)^. Thus, if we do not restrict the range of Q^, we 
obtain the same sufficient condition for stability as (|6-9|) , but for the restricted range 

< 8n(n — l)M2/(3n — 2)^, we obtain the stronger sufficient condition ()6-ll() . 
which coincides with the condition obtained in Paper II for the case Q = 0. 

Similarly, for K = —1 and A < 0, rewriting Vr as 
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n+1 





K=1 




A,<0 


V ^ H max 
\ (^=^cl) 










X>0 






^ 








1 (n+1)2 
4n 





K=0,-1 




1 



Fig. 2. Ranges of the value xh- 



we obtain a sufficient condition for stability stronger than (|6-9() . 

XL + 3n-2 = k^ + n-2>0, 
from Vt > if we restrict the range of A to 



(6-13) 

This condition is sufficient to guarantee the stability of a maximally symmetric black 
hole with K = —1 for n > 2. However, if we extend the range of A to the whole 
allowed range, i.e., that satisfying A > Ac-, (|(i-9j) is the strongest condition that can 
be obtained only from Vr > 0- 



6.2. Vector perturbation 

For V = Vv± in (11^ and 



5 = ^, 
2r ' 



we obtain 



my + 



l)M ± A 



„n— 1 



(6-15) 

(6-16) 
(6-17) 



my = ky — [n — 1)K. 

Hence, from ky > [n — 1)K, we find that Vv+ is always positive, and static charged 
black holes are stable with respect to the electromagnetic mode of the vector per- 
turbation. In contrast, from 



Vv- 
h := 



f 



myh 



r2 (n2 - 1)M + Z\' 



, p , 2n(n - l)g2 



(6-18) 
(6-19) 
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Vs. 









..---^n=5 




n=4 




. n=3 




n=2 


1 // " 


2 


^ r/rH 



K=1, ;i<0, Q^/M^=0.1, 1=2 



Fig. 3. Examples of Vs+ for A" = 1 and A < 0. 

it is seen that Vy- may become negative: Because /i is a monotonically increasing 
function of r, Vy- is positive if and only if hiju) > 0. 

6.2.1. The case A > 

In this case, the background spacetime contains a regular black hole only for 
K = 1, and the static region outside the black hole is given hy rn < r < rc (< +cxd). 
In this region, from (|A-29b|) . we have 2M/r"'"^ < xh < 2;//,max < (n + l)M^/(n(5^). 
Hence, h > 0, and the black hole is stable in this case. 

6.2.2. The case A < 

In this case, as shown in Appendix the spacetime contains a regular black 
hole if A > Ac and 2M/r"^^ < xh < 3;_f/,max- Hence, from (|A-29b|) . we obtain the 
relation 

h > V(n2 - 1)2M2 + 2n(n - l)mvQ'^ - \l {n^ - ifW^ - AKn{n - l^Q^. (6-20) 

For K = 0,1, h > follows from this. Hence, the black hole is stable. In contrast, 
for K = —1, the right-hand side of this inequality becomes negative for ky < n — 1. 
Hence, Vy- can become negative near the horizon if A is sufficiently close to Ac-, 
provided that the spectrum of ky extends to A:^ < n — 1. 

6.3. Scalar perturbation 

By applying the S'-deformation to Vs+ with 

S=f^; h^ = r-/'~^H^, (6-21) 
n+ ar 

we obtain ^ 

^s+ = ^fj^[{n-2){n + l)5x + 2]. (6-22) 

Since this is positive definite, the electromagnetic mode ^4. is always stable for any 
values of K, M, Q and A, provided that the spacetime contains a regular black hole, 



32 




Fig. 4. Examples of Vs-. 



although Vs-\^ has a negative region near the horizon when A < and Q'^/M'^ is smah 
(see Fig. EJ. 

Using a similar transformation, we can also prove the stability of the gravita- 
tional mode ^_ for some special cases. For example, the ^-deformation of V^- with 

h- ar 

leads to 

^S- = iTTTT- - (n + l)(n - 2)(1 + m5)x] . (6-24) 
2r^/i_ 

For n = 2, this is positive definite for m > 0. When = 1, A > and n = 3 or 
when A > 0, Q = and the horizon is 5^, from m > n + 2 > 2) and the behaviour 
of Xh (see Fig. we can show that V^. > 0. Hence, in these special cases, the 
black hole is stable with respect to any type of perturbation. 

However, for the other cases, Vs- is not positive definite for generic values of 
the parameters. The S'-deformation used to prove the stability of neutral black holes 
in Paper II is not effective either. This is because Vs- has a negative region around 
the horizon for the extremal and near extremal cases, as shown in Fig. [IJ and the 
S'-deformation cannot remove this negative region if 5 is a regular function at the 
horizon. Hence, determination of the stability for these generic cases with n > 3 is 
left as an open problem. 

§7. Summary and Discussion 

In the present paper, we have extended the formulation for perturbations of a 
generalised static black hole given in Paper I to an Einstein-Maxwell system in a 
generalised static spacetime with a static electric field, and we have shown that the 
perturbation equations for vector and scalar perturbations can be reduced to two 
decoupled second-order ODEs for the gravitational mode and the electromagnetic 



33 



Table I. Stabilities of generalised static black holes. In this table, "d" represents the spacetime 
dimension, n + 2. The results for tensor perturbations apply only for maximally symmetric 
black holes, while those for vector and scalar perturbations are valid for black holes with generic 
Einstein horizons, except in the case with K = 1,Q = 0, X > and d — 6. 







Tensor 


Vector 


Scalar 






Q = 




Q = 




= 


Q/0 


K = 1 


A = 


OK 


OK 


OK 


OK 


OK 


d = 4, 5 OK 
d> 6 ? 




A > 


OK 


OK 


OK 


OK 


rf < 6 OK 
d > 7 ? 


d = 4, 5 OK 
d> 6 ? 




A < 


OK 


OK 


OK 


OK 


d = 4 OK 
d > 5 ? 


d = 4 OK 

d > 5 ? 


K ^0 


A < 


OK 


OK 


OK 


OK 


d = 4 OK 
d > 5 ? 


d = 4 OK 
d > 5 ? 


K = -1 


A < 


OK 


? 


OK 


? 


d = 4 OK 
d > 5 ? 


d = 4 OK 
d > 5 ? 



mode, irrespective of the value of the cosmological constant and the curvature of the 
horizon. In particular, we have found that the coupling between the perturbations of 
the metric and the electromagnetic field produces significant modifications of the ef- 
fective potentials for the gravitational mode and the electromagnetic mode when the 
black hole is charged. Our formulation also provides an extension of the correspond- 
ing formulation for the Reissner-Nordstrom black hole in four dimension^^''^''^ 
to asymptotically de Sitter and anti-de Sitter cases. 

With the help of this formulation and the method used in Paper II, we have 
analysed the stability of generalised static black holes with charge. The results are 
summarised in Tabled As shown there, maximally symmetric black holes are stable 
with respect to tensor and vector perturbations over almost the entire parameter 
range; the exceptional case corresponds to a rather exotic black hole, whose horizon 
is a hyperbolic space. In contrast, for a scalar perturbation, we were not able to 
prove even the stability of asymptotically flat black holes with charge in generic 
dimensions, due to the existence of a negative region in the effective potential around 
the horizon in the extremal and near extremal cases, in contrast to the neutral case. 
Whether this negative ditch produces an unstable mode or not is uncertain. Hence, 
the stability of asymptotically flat and asymptotically de Sitter black holes for d > 6 
and of asymptotically anti-de Sitter black holes for d > 5 are left as open problems. 
In connection to this, it should be noted that the existence of a negative region 
in the effective potential near the horizon may have a significant influence on the 
frequencies of the quasi-normal modes and the greybody factor for the Hawking 
process, even if these black holes are found to be stable. 

In the present paper, we have also given explicit expressions for the source terms 
in the master equations. As mentioned in the introduction, this information will be 
necessary in the estimation of gravitational and electromagnetic emission from black 
holes in higher dimensions. In addition to this practical application, we also expect 
that these master equations with source terms can be used in the analysis of static 
singular perturbations of black holes associated with some singular source such as 
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a string or a membrane. For example, if one treats the C-metric as a perturbation 
of a spherically symmetric solution, it is found that this perturbation obeys an 
equation with a string source. Hence, it is expected that one can obtain some 
information concerning higher-dimensional analogue of the C-metric by studying 
singular solutions to the master equation with a singular source. 
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Appendix A 

Parameter range for the existence of a regular black hole 

In this section, we determine the parameter range in which the background 
metric ()2-18() contains a regular black hole. Here, by a regular black hole, we mean a 
degenerate Killing horizon or bifurcating Killing horizons at r = rjj that separate(s) 
a regular region with /(r) > and a singular region containing the singularity at 
r = 0. 

First, we consider parameter values satisfying A > and K < 0. In this range, 
if Q = 0, there exists no horizon, because / is negative everywhere. If Q ^ 0, from 

(r"-V)' ={n- l)i^r"-2 - (n + l)Ar'" - ~ ^^^^ < 0, (A-1) 



/ is a monotonically decreasing function of r. Further, 



"'1 <-2M, r^+oo • ^^^^ 

Hence, the spacetime has no region with / > that is separated from the singular 
region by Killing horizons. Therefore, the spacetime contains regular black holes 
only for ii' = 1 or A < 0. 

A.l. The case Q = 

In this case, since / —oo as r ^ +0, the spacetime contains a regular black 
hole if there is a region in which / > 0. From 

./ , 2(n-l)M , 

/ = -2Ar+ \J , (A-3) 

it is seen that the result depends on A. 
A.1.1. A < or = 1 and A = 

Because / is monotonically increasing and becomes positive as r ^ +oo, / has 
a single zero at r = rn, and / > for r > rn- Thus, there is a regular black hole. 
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A. 1.2. K = 1 and A > 

Because / has a single maximum at r = a, the spacetime contains a regular 
black hole if and only if /(a) > 0. Then, since f'{a) = is equivalent to 

A='"-';^. (A.4) 

this condition can be written as 

a"-^ > (n + 1)M. (A-5) 

In terms of A, it is expressed as 

-. (A-6) 

f has two zeros, ai r = rn and r = Vc {th < fc), and / > for rn < r < Vc, 

A.2. The case Q^O 

When Q ^ 0, since / +oo as r ^ +0, there must exist a point r = a such 
that /'(a) = and /(a) < in order for the spacetime to contain a regular black 
hole. From 

/' = -2Ar + 2(n-l)(^ (A-T) 

A is expressed in terms of a as 

, , (n - 1) / M Q2 x 

From this, /(a) can be written as 

^/ ^ (n + 1)M , raQ2 

Hence, the condition /(a) < requires 

D:= {n+ ifNP - AKnQ^ > 0, (A-10) 
and under this condition, /(a) < if a satisfies 

(n + l)M-/D\ 1 ^{n + l)M + ^/D 

max 0, -^^ —-7^ < ^ < —-7^ . (A-11) 

\^ ' 2nQ2 j - a*^-! - 2nQ^ ^ ' 

Note that g{a) has a maximum at a = a^, where 



and it is monotonic everywhere except at this point. Further, ^(+0) = — oo and 
5(+oo) = 0. Hence, A < for a < oq, where 

= (A- 13) 

and A > for a > qq. 
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A. 2.1. K = 1 and A = 

In this case, the spacetime contains a regular black hole if and only if > Q^, 
and the horizon is at r"~^ = r^~^ = M + ^/M^ - Q^. 

K.2.2. if = 1 and A > 

Because the sign of /'(r) is the same as the sign of g{r) — A, the condition that 
there is a point r = a such that /'(a) = is equivalent to the relation 

A < A„,a^ := g{a^) = ^ ^ ^ ( ) ( ^ ) —. (A-14) 

M"-i 



Further, from the condition D > 0, we have 

Under these conditions, /'(a) = has in general two solutions, a = ai,a2 (ao < 
O'l < CL-m < «2), and the spacetime contains a regular black hole if and only if 
/(ai) < < /(as). 

First, from /(as) > 0, we obtain 

In terms of A, this can be expressed as 

, , , , ^n+in-1 {n-l)M + y/D . ^ ^. 

A < \c2 ■■= g{ac2) = 2"-i ^ '- — (A-17) 

[{n + 1)M + VD]"^ 

Next, because /(a) becomes minimal at a = a^, from /(oq) = 1 — M'^ /Q'^, we have 
/(oi) < /(oo) < if Q2 < Also, for < Q"^ < {n + lf/{An)M^, f{ai) < is 
equivalent to 

a"- < =: a"-, (A-IS) 

or, in terms of A, to 

A>Aci :=5(aci) = 2"-i ^ (A-19) 

^ [(n + l)M-/D]~ 

Here, if we vary with M fixed, we have 

-r-{K2 - Aei) = ^ K+2 - ^ < 0. (A-20) 

Hence, Ac2 — Ad is a monotonically increasing function of Q for fixed M, and Ad = 
Ac2 at = (n + l)^M^/(4n). From these results, it follows that Ad < \c2 for 
Q2 < (n + l)^/(4n)M2. Further, Ad vanishes at Q"^ = M'^. Therefore, for A > and 
K = 1, the spacetime contains a regular black hole if and only if 

max (0, Ad) < A < A^a : < ^ ^ (A-21) 

47T, 
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Table II. The parameter range in which the spacetime contains a regular black hole. 





K = 1 


K = 


K = -1 


A = 


cr < 






A > 


For < M'^ 

A < \c2 

For <Q^ <{n+ l)^M'^/{4n) 
Aci < A < Ac2 






A < 


Q'^ < M'-' and Aci < A 


Aco < A 


Ac- < A 



A.2.3. A<0 

In this case, /'(a) = always has a single solution, and / — +oo as r — +0, +oo. 
Hence, the spacetime contains a regular black hole if and only if /(a) < 0. Because 
a < ao < Qm, this condition is equivalent to 

J_ < J_ = jn + m + VD 

or in terms of A, 

This, together with D > 0, leads to the following conditions: 

^ = 0: A >A.o :=-—(— ) (^) (A.24a) 

i. = -l:A>A._:=-2A!L^-^^(=-il^i^, (A.24b) 

" [V^-(n + l)M]^ 

. > X VD-(n-l)M x 
i^ = l: A>Aei:=-2n-i ^ (A-24c) 

" [{n + l)M-VD]^ 

Finally, we determine the range of xh = 2M/r^~^. In general, A is expressed in 
terms of rn as 



Prom this, we have 



dVH 

First, for iiT = 1, this can be written in terms of and as 
dA _ 2nQ^ / ]_\ / 



(A-27) 
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From this and the relation ac2 < a2 < rn < ai < ad, it follows that A is a mono- 
tonically increasing function of rn for fixed M and Q. Hence, from the constraint 
Aci < A < Ac2, we obtain 

,max) 

(A-28) 

where XH,m\n and XH,ma.^ are the values of xh for A = Ac2 and A = Ad, respectively: 



4M 



XH,i 



XH,i 



(n + l)M + ^/D' 
(n + 1)M2 + My/D 



Next, for K = OT K 
dX 2nQ2 



-1, dX/drn can be written in terms of ac as 
1 1 \ 



dru 



' H 



1 _^ - (n + 1)M 



(A-29a) 
(A-29b) 

(A-30) 



' H 



' H 



2nQ2 J ■ 

From rj{<a< a^, this is non-negative. Hence, from Ac < A < 0, we obtain 

< XH < XH,raa.x- (A-31) 



The allowed ranges of xh given in (|A-28|) and ()A-3H) in the {xh,Q'^/M'^) plane are 
displayed in Fig. [2 

Appendix B 



Expressions for Ej and El 

We have the following expressions for Ef and El- 

2Ei = -fdliX + St) + - ^) dr{X + St) 



+ 



P + {n-2)K (n-2)(n-3)M 



+ 



-3(n- 1)A + 



(2n-3)Q 



21 



r,2n 



{X + St) 



-fd',Y-{^ + ^]drY 



+ 



2K (n-2)(n + 3)M (n + l)(2n - 3)Q2n 
h 3A H 



r.2n 



(B-1) 



1 



2El = jd^X + ^drX 



2 



+ 



(n-l)(n-2)M n(n-l)Q2n 
— — 3A H — 



„2n 
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+ 



(n - l)(n - 2)M ^ _ 3(n- l)(n - 2)Q^^ 



r.2n 



Y 



2 



2 



+ -dtdrZ +{—-^]dtZ 



+ 



(n-l)P 2(n-l)(n-2)M 



+ 



-3(n- 1)A 



2(n- l)(n- 3)Q- 



21 



Appendix C 

Expressions for coefficient functions 



We have the following: 

Pg^ = [-4n^z + 2n^(n + l)x - 4n(n - 2)m] y 

+ {2n^(n - l)x + An{n - 2)m + 4n^(n - 2)i^} z 
-n^(n^ - l)x2 + {-4n(n - 2)m + 2n^(n + l)K] . 



(B-2) 



(C-1) 



+4m^ + An^mK, 
Ps2 = \pn^z — n^(n + l)(n + 2)x + 2n{n — 4)m] y 

-2n^z^ + {n^(3n^ - n + 2)x - 4ri(ra - 2)m - 6n^(ra - 1)/^} z 
-n^{n + + {n(3n - 7)m + n^(n^ - \)K] x 

-2m2 - 2n{n - l)mK, (C-2) 

P53 = -2n^(3n - 2)z + n^{n + l)x - 2(n - 2)m. (C-3) 



= [4n^{n^ - l){n^z + m)a; - 8n^{n - l)(3n - 2)z^ 
-32n^(n - 1) Vz - 8(n - l)(n - 2)m^] y 
+n^(n + l)-lx''^ + 2n(n + 1) {-n2(3n - l)(n + 2)z 

+2(n^ + n + 2)m - n(n - 2)(n + 
+ [4n''(13n - T)^^ - 4n^ {(971^ - 8n + 15)m + n(n + l)(n^ - 7n + 8)K} z 

-An{n - ll)m^ - 4n^(n + l)(n - 5)mK] x 
-8n^(3n - 2)z^ + Sn^ {4(n2 _ n + l)m + n2(3n - 4){n - 2)K] z^ 
+ {(16n^ - 72n + \%)rr? + 16n^(2n - 3)mii:} z 

+8n^m^K + IGm^, (C-4) 
Pxi = [4n^(n - 1)2; + 4(n - l)m] y 
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+n(n + l)lx^ + {-8n^z + 2(3n - l)m - 2n(n + 1)K} 
+Av?z^ - {4(2n - l)m + 4n^(n - 2)K] z - inmK, 



X 



(C-5) 



Pyo = [2n^{n + ifx"^ - 4n^{n + 1) {n^(5n - 3)z + {n - 3)m} x 

+8n'^(2n - l)(3n - 2)z^ + 16n^(n^ - 4n + 2)mz - 8(n - 2)m^] y 
+n^(n - l)(n + 

+2n(n^ - 1) {-2n^(3n - l)z + 4m - n(n - 2)(n + l)K} 

+ [20n^(2n - l)(n - 1)^^ + 4n2(n - 1) {(n - 13)m + n(5n - 8)(n + l)i^} z 

+12ra(n - l)m^ + 4n^(n^ - l)mEr] x 
-8n^(3n - 2)(n - l)z^ + 16n2(n - 1) {2m - n^(3n - 4)ii:} 
-8(n - l)m { (5n - 2)m + 2n^ii:} (C-6) 
Pvi = [2n^(n + l)x - 4n^(2n - l)z + 4m] y 

+n(n^ - - 2(n - 1) {Sri^z + m + n(n + l)if} x 

+An'^{n - \)z^ + 4(n - l)(m + 2n^ii:)z, (C-7) 
Pz = [-n^(n + l)x + 2n^(3n - 2)z + 2(n - 2)m] y 

+n{n + + {n^(3n - l)z + (4n - 2)m + n{n + l)(n - 2)K] x 

-2r?{n - 2)^2 - {(6n - 4)m + 2n^(3n - A)K] z - 2nmK. (C-8) 

Note that these functions satisfy the following relations: 



Pxo + Pyo + ^nHPz = 16H^, 
Pxi + Pyi = -4nrfH. 



(C-9a) 
(C-9b) 



PxA = [4n^(n - l)z + 4n{n - l)m] y 

+n^(n + ifx"^ + 4n {-2n^z + (n - l)m - n(n + 
+4n^z'^ - 4n { (n - l)m + n^(n - 3)K} z 
-Arri^ - 4n(n + l)mK, 



X 



(C-10) 



P^2 = [4(n - l)n^z + 4(n - l)m] y 

+n(n + ifx'^ + {-8n^z + (6n - 2)m - 2n(n + 1)K} 
+4n^z^ - {4(2n - l)m + 4n^(ra - 2)^:} z - AnmK, 



X 



(C-ll) 



P^3 = (n + l)x + 2n(n - 2)z + 2m, 



(C-12) 



Pya = [2n^(n + - 4n^(2n - l)z + 4nm] y 

W(n^ - + {-6n^(n - 1)2; + 4nm - 2n^(n + l)(n - 2)K} x 
+4n^{n - \)z^ + 4n {-m + r?{2n - 3)K] z + 4nmK + 4m^, (C-13) 
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Py2 = [2r?{n + l)x - 4n^(2n - l)z + 4m] y 

+n(n^ — l)x^ — 2(n — 1) {Sn^z + m + n{n + 1)-?^} x 

+4n^(n - 1)^2 + 4(n - l)(m + 2n2K)z, (C-14) 

Pyg = (n - 1) [(n + l)x - 4nz] . (C-15) 
The following relations hold for these functions: 

PxA + PYA = -4.n^fH, Px2 + PY2 = -^nfH, Px3 + Py3 = 2H. (C-16) 

Appendix D 

Gauge-invariant treatment of the exceptional modes 



In the present paper, we imposed the Einstein equation + 2(n — 2)F = 
as the gauge condition for the exceptional modes of the scalar perturbation with 
k"^ = n for K = 1. As discussed in Paper I, this gauge condition does not fix the 
gauge freedom completely, and it leaves the residual gauge freedom represented by 
5z* = L§* satisfying 

. (r-DL) + = 0. (D-l) 

Hence, the master variables introduced in ^ contain unphysical degrees of freedom 
for the exceptional modes, although the master equations themselves are gauge- 
invariant. In this appendix, we express the master equations in terms of genuinely 
gauge-invariant variables. 

We define all quantities for an exceptional mode corresponding to the gauge- 
invariant variables for a generic mode introduced in the text by the same expressions 
with = 0. Then, for the general gauge transformation ()2-31|) . Xa transforms as 

6Xa with 

6Xa = -r^Da{L/k)+Ta. (D-2) 
From this, the following transformation laws of F and Fj^ are obtained: 

F^F + 6F: 6F = -rDr ■ D{L/k) - L/k, (D-3a) 

Fab ^ Fab + SFab : SF^b = -Da [r^D.j^ - l^r^Daj^ . (D-3b) 

In particular, we have 

6{F^ - 2nF) = -2 

Hence, from (|2-32|) and (|5-1()|) . £ and £a transform as 



V r J ^ 



(D-4) 



5E = -D. {r^E,D{L/k)) = -qD • (^^^^) , (D-5a) 
S£a = rEoDaL = ^DaL. (D-5b) 
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From this and the definition of A, (|5-19|) . we find that A transforms as 

L 



6A = q 



k' 



(D-6) 



Like these variables, the matter variables Uab^^a and Ul are gauge dependent. 

However, Sa,S^ and Sl are gauge-invariant, because they represent perturbations 

of quantities whose background values vanish, like Ja- 
in order to proceed further, we have to treat the black-hole-type background 

case and the Nariai-type background case separately. 

D.l. Black hole background 

In the black hole background (|2-18)) , the gauge transformation of is written 



SFI 



2wVL 



Tl 



■' k V 2 k 



-2r'fiL"/k)-r{Af + rf'){L'/k). 



In particular, we have 



6{2F 



2HL 

n k 



k 



Hence, the master variable ^ transforms as 

For k"^ = n and IT = 1, ^ = Af and <P± are written 



{n + l)MQ / <P 2A\ 

j.n/2-l \^^n/2 ^ q ) ' 

2{n + 1)M 1 ^ ^ 



„n— 1 



r.nl2-\ 



(D-7a) 

(D-7b) 
(D-7c) 

(D-8) 
(D-9) 

(D-lOa) 
(D-lOb) 



From this, we find that ^+ is gauge- invariant. In contrast, as is seen from the relation 

^= (D-ll) 

2(n + l)M (n+l)M ^ ' 

is not gauge- invariant. Nevertheless, the master equation for is gauge- 
invariant. This becomes evident if we rewrite this equation in terms of the gauge- 
invariant combinations 



F*- J- 



nrf 



f 



Ujjr 

,n/2 



(D-12a) 
(D-12b) 
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pr _ pr 



1 



+ 



j.n/2 

n{n — 2)f — nrf 



F := F 



1 



2rn/2 



rf<P' + 1 



n/ 
2 



(D-12c) 
(D-12d) 



Then, taking account of the fact that the master equation was derived under the 
gauge condition F^ + 2(n — 2)F = 0, we obtain 

1 



F^ + 2{n-2)F 



2(n+ 1)MW2-2/ 



Here, for k'^ = n and K = 1, we have 



+ 2nVdSs+ + 2nVdiVs+ - Vs^)<P^ 

(D.13) 



- K- = ^rr^^ [ {4 - 2n\n + l)5x} y 

+6{n'^ - + {2n(n + l)(n - 2)J - 2n + 2}2;]. 



(D-M) 

Since this equation is written only in terms of gauge-invariant variables, it is valid 
in any gauge. Thus, the master equation for ^_ gives an algebraic relation among 
the gauge- invariant variables F and F^. 
The definition of <P yields another relation, 

2iujrF + F[ = 0. (D-15) 

Further, the expressions H5-45() for X and Z provide two more relations: 



„n-2 



{H - 2F) 



v?{n - l)Q ( nPx 



-lOJ- 



2{n + l)Mr"/2+i V 
n2(n- 



If''- 



+ Z. 



sU=Oi 



+ X,U=o,(D-16a) 
(D-16b) 



2n{n-4)z. (D-17) 



(n-hl)MW2i?' 

where 

Px^ = [{n + l)(n - 2)x + 2n^z\ y^(n^ \){2n - l)^^ 
-{(3n2 3n - 2)z (n -h l)(n 2)}x 2n^z-^ ■ 

The four equations dTTTIl) . (ITTTKl) and (ITTT6|) can be solved to yield expressions 
for F and -F^ in terms of ^+ and the gauge- invariant matter source terms 5a, S*^ 
and Jfl. Therefore, the only dynamical gauge- invariant variable for the exceptional 
modes is Note that we can impose the gauge condition ^ = 0, and for this 
gauge, F^ = F^ and F = F hold and is proportional to A. 

D.2. Nariai-type background 

The case of the Nariai-type background can be treated in almost the same way. 
First, the gauge transformations of F^ and F are written 



6Fl 



-2a^ 



y ^ '2T 



(D-18a) 
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-SF,^ = 2ia'c.{f^-^^], 



k 



dF = - 



2 k 

I Tl 



1 ( L" f L 
L 



k' 

For m = and K = 1, a and f^ have the simple expressions 

n — 1 n{n — 1)Q'^ 



a = 



Hence, <P± are written 



n - 1 _^ n(n- 1)Q2 



2n 



2(n-l)Q 



2(n- 1) 



a 



— 

Q2n-2 



(D-18b) 
(D-18c) 
(D-18d) 

(D-19a) 
(D-19b) 

(D-20a) 
(D-20b) 



Prom this, we find that is gauge-invariant, while is not gauge-invariant, as 
seen from the relation 



nQ 



$4 



(D-21) 



as in the black hole background case. 

We can construct the following gauge-invariant qTiantities from F^ and F: 



^F + ^F' 
f 2 



F^ = Fl - 2a^ 



Ff = Ff + 2iua^ l^fF'-^F 
pP = FP - 2c? ijF" + ^F'^ . 



Under the gauge condition F^ -\- 2(n — 2)F = 0, we have 



n 



(D-22a) 
(D-22b) 
(D-22c) 

(D-23) 



From this, it follows that the master equation for can be expressed in terms of 
the gauge invariant variables as 



F^ 



m7 



nQ 

,2n-4^+" 



Further, the equations for X and Z can be rewritten as 



nQ 



-SP 



na 



n— 4 * 



(D-24) 

(D-25a) 
(D-25b) 
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The corresponding expression for Fp can be obtained from those for and F^. 
Under the gauge condition F = 0, F^ coincides with F^, and ^+ becomes a constant 
multiple of A. 
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